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Abstract 

Charm production in deep inelastic scattering has been measured with the ZEUS 
detector at HERA using an integrated luminosity of 82 pb _1 . Charm has been 
tagged by reconstructing D* + , D°, D + and (+ c.c.) charm mesons. The 
charm hadrons were measured in the kinematic range px(D* + , D°, D + ) > 3 GeV, 
Pt{D+) > 2 GeV and \rj(D)\ < 1.6 for 1.5 < Q 2 < 1000 GeV 2 and 0.02 < y < 0.7. 
The production cross sections were used to extract charm fragmentation ratios 
and the fraction of c quarks hadronising into a particular charm meson in the 
kinematic range considered. The cross sections were compared to the predictions 
of next-to-leading-order QCD, and extrapolated to the full kinematic region in 
Pt{D) and r](D) in order to determine the open-charm contribution, F 2 cc (x, Q 2 ), 
to the proton structure function F 2 . 
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1 Introduction 



Heavy-quark production in ep interactions in the deep inelastic scattering (DIS) regime is 
dominated by the interaction between the exchanged virtual photon and a gluon within 
the proton, the so-called Boson Gluon Fusion (BGF) mechanism. Heavy-quark production 
provides a twofold test of perturbative quantum chromodinamics (pQCD): a study of the 
BGF process and the higher order corrections to it, and an independent check of the 
validity of the gluon density in the proton extracted from the inclusive DIS data. Of the 
two heavy quarks whose production is accessible by HERA, c and b, the latter is strongly 
suppressed due to its smaller electric charge and larger mass. This paper reports a study 
of c-quark production. 

A charm quark in the final state is identified by the presence of a corresponding charmed 
hadron. This paper studies the production of the pseudo-scalar mesons D°, D + , and 
the vector meson D* + from the decays D° — > K~ir + , D + — > K~tt + tc + , — + <f)ir + — ► 
K + K~ir + and D* + — > D°tt + — > K^-k + -k + (the charge conjugated modes are implied 
throughout this paper). Since a D hadron is measured and not the c quark itself, any 
comparison with pQCD requires a modelling of the c — > D fragmentation. A consequence 
of the QCD factorisation theorem [1] is that the "hard" (pQCD governed) c-production 
mechanism is independent of the "soft" fragmentation process. Measurements of D- 
hadron cross sections provide therefore information about both c-quark production and 
its fragmentation. 

This paper presents a complete study of D-meson production in DIS at HERA: measure- 
ments of c-quark fragmentation ratios and fractions with unprecedent precision, D-meson 
differential cross sections and the charm contribution, F 2 CC , to the proton structure func- 
tion F 2 . It addresses the universality of fragmentation and tests the predictions of pQCD 
for charm production. The data sample used was taken by the ZEUS detector during the 
years 1998 - 2000. The fragmentation measurements follow closely those reported recently 
by ZEUS in the photoproduction regime [2]. Using a variety of D mesons, the pQCD anal- 
ysis complements the study done with D* + in the same data sample [3]. Measurements 
of D* cross sections are only used in this paper for the extraction of the fragmentation 
parameters. 

Similar measurements of the properties of c-quark fragmentation in DIS have also been 
performed by the HI collaboration [4]. Other previous measurements of charm production 
in DIS with pQCD analyses used the D* + meson [3,5-9] or inclusive lifetime tags [10]. 
There are also several measurements of charm photoproduction [9,11-15]. 
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2 Experimental set-up 



The analysis was performed with data taken from 1998 to 2000, when HERA collided 
electrons or positrons with energy E e = 27.5 GeV on protons of energy E p = 920 GeV. 
The results are based on e~p and e + p samples corresponding to integrated luminosities 
of 16.7 ± 0.3 pb~ x and 65.1 ±1.5 pb" 1 , respectivelyJl] 

A detailed description of the ZEUS detector can be found elsewhere [16]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles are tracked in the central tracking detector (CTD) [17], which oper- 
ates in a magnetic field of 1.43 T provided by a thin superconducting coil. The CTD 
consists of 72 cylindrical drift chamber layers, organised in 9 superlayers covering the 
polar-angl^l region 15° < 6 < 164°. The transverse-momentum resolution for full-length 
tracks is a(p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-scintillator calorimeter (CAL) [18] consists of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part 
is subdivided transversely into towers and longitudinally into one electromagnetic sec- 
tion and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections. The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a(E)/E = 0.18 j\[E for electrons and 
a(E)/E = 0.35/ \[E for hadrons, with E in GeV. 

The position of the scattered electron was determined by combining information from 
the CAL, the small-angle rear tracking detector (SRTD) [19] and the hadron-electron 
separator (HES) [20]. 

The luminosity was measured from the rate of the bremsstrahlung process ep —>■ ejp, 
where the photon was measured in a lead-scintillator calorimeter [21] placed in the HERA 
tunnel at Z = —107 m. 

3 Theoretical predictions 

The next-to-leading order (NLO) QCD predictions for cc cross sections were obtained us- 
ing the HVQDIS program [22] based on the so-called fixed-flavour-number scheme (FFNS). 
In this scheme, only light partons (u, d, s, g) are included in the proton parton density 

1 Hereafter, both electrons and positrons are referred to as electrons, unless explicitly stated otherwise. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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functions (PDFs) which obey the DGLAP equations [23] , and the cc pair is produced via 
the BGF mechanism [24] with NLO corrections [25]. This calculation is expected to be 
valid [26] in the kinematic range of this measurement, 1.5 < Q 2 < 1000 GeV 2 , where Q 2 
is the negative of the four-momentum transfer squared, hereafter referred to as photon 
virtuality. 

The following inputs have been used to obtain the predictions for .D-meson production at 
NLO using the program HVQDIS. The FFNS variant of the ZEUS-S NLO QCD fit [27] 
to structure-function data was used as the parameterisation of the proton PDFs. In this 

(3) 

fit the three-flavour QCD scale Aqcd was set to Aq^ d = 0.363 GeV and the mass of 
the charm quark was set to 1.35 GeV; the same mass and Aq^ d were therefore used 
in the HVQDIS calculation. The renormalisation and factorisation scales were set to 
\x = \/Q 2 + Am 2 for charm production both in the fit and in the HVQDIS calculation. 
The charm fragmentation to the particular D meson was carried out using the Peterson 
function [28]. The values used for the hadronisation fractions to D mesons, /(c — > D), 
were those measured in this paper, and the Peterson parameter, e, was set to 0.035 [29]. 
The effect of J/ip production was found to be negligible [30,31]. 



4 Kinematic reconstruction and event selection 

The kinematic variables Q 2 , the Bjorken scaling variable, x (in the quark-parton model 
x can be interpreted as the fraction of proton momentum carried by the struck quark), 
and the fraction of the electron energy transferred to the proton in the rest frame of the 
proton, y, can be reconstructed using a variety of methods, whose accuracy depends on 
the variable of interest and its range: 

• for the electron method (specified with the subscript e), the measured energy and 
angle of the scattered electron are used; 

• the double angle (DA) method [32] relies on the angles of the scattered electron and 
of the hadronic system; 

• the Jacquet-Blondel ( JB) method [33] is based entirely on measurements of the hadronic 
system; 

• the E-method [34] uses both the scattered-electron energy and angle, and measure- 
ments of the hadronic system. 

The reconstruction of Q 2 and x was performed using the S-method, since it has better 
resolution at low Q 2 than the DA method. At high Q 2 , the S-method and the DA method 
are similar, and both have better resolution than the electron method. The DA method 
was used as a systematic check. 
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A three-level trigger system was used to select events online [16,35]. At the third level, 
events having at least a reconstructed D* + ,D°,D + (only e + p sample), D£ or Ajt can- 
didate, as well as a scattered-electron candidate, were kept for further analysis. The 
efficiency of the online reconstruction for any of the above hadrons, determined relative 
to an inclusive DIS trigger, was generally above 95%. 

The events were selected offline [5,36] using the following cuts: 

• the scattered electron was identified using a neural-network procedure [37,38]. Its 
energy, E e >, was required to be larger than 10 GeV; 

• 2/ e <0.95 (where y e is y reconstructed with the electron method) and ?/jb > 0.02 (where 
?/jb is y reconstructed with the JB method). The former condition removes events 
where fake electrons are found in the FCAL and the latter rejects events where the 
hadronic system cannot be measured precisely, in order to reconstruct the kinematic 
variables; 

• 40 < 5 < 65 GeV, where 5 = J^-Ei(l — cos^j) and Ei and 0, are the energy and 
the polar angle of the i th energy-flow object (EFO) [39] reconstructed from charged 
tracks, as measured in the CTD, and energy clusters measured in the CAL. The sum 
i runs over all EFOs; 

• a primary vertex position determined from the tracks fitted to the vertex in the range 

| ^vertex | < 50 Cm; 

• the impact point (X, Y) of the scattered electron on the RCAL was required to lie 
outside the region 26 x 14 cm 2 centred on X = Y = 0. 

The angle of the scattered electron was determined using either its impact position on the 
CAL inner face or a reconstructed track in the CTD. When available, SRTD and HES 
were also used. The energy of the scattered electron was corrected for non-uniformity 
effects caused by cell and module boundaries. 

The selected kinematic region was 1.5 < Q 2 < 1000 GeV 2 and 0.02 < y < 0.7. 

5 Reconstruction of charm hadrons 

The production of D* + , D°, D + , and Df charm mesons was measured in the range of 
transverse momentum Pt{D) > 3 GeV and pseudorapidity |?7(-D)| < 1.6. For the D+, 
the pt{D^) requirement was relaxed to px{Df) > 2 GeV, as the constraint provided by 
the K + K~ pair coming from a <fi meson kept the combinatorial background at acceptable 
levels. The reconstruction of the baryon was attempted using the decay A+ — > K~pn + . 
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The signal achieved had a statistical significance of around three standard deviations, and 
therefore it was not used. 

The charm mesons were reconstructed using tracks measured in the CTD and assigned 
to the reconstructed event vertex. To ensure good momentum resolution, each track 
was required to reach at least the third superlayer of the CTD. Further background 
reduction was achieved by imposing cuts on the transverse momenta and decay angles 
of the charm-hadron decay products. The cut values were optimised using Monte Carlo 
(MC) simulation to enhance signal over background ratios while keeping acceptances high. 

The cross sections of the D° and D* + mesons and related quantities involved in the 
measurements of fragmentation properties can be obtained from the combination of three 
independent samples [2]: those of D° candidates with and without a "AM" tag and that 
of "additional" D* + candidates. The samples are described below. The rationale for this 
division [2] will become apparent in Sections M and M 

5.1 Reconstruction of D° mesons 

The D° mesons were reconstructed from the decay channel D° — > K~ir + . In each event, 
D° candidates were formed from pairs of tracks with opposite charges and px > 0.8 GeV. 
The nominal kaon and pion masses were assumed in turn for each track. To reduce the 
combinatorial background, a further cut was applied in the angle, 9*(K), between the 
kaon in the D-candidate rest frame and the D-candidate line of flight in the laboratory 
frame, | cos 6*(K)\ < 0.85. The D° candidates were separated into two groups. The AM 
tag group consists of D° candidates that, when combined with a third track that could be 
a "soft" pion (vr s ) in a D* + -> D°ir+ decay, have AM = M(Ktxtx s ) - M (Kn) in the range 
0.143 < AM < 0.148 GeV. The soft pion was required to have px > 0.2 GeV and charge 
opposite to that of the kaon. This cut was raised to px > 0.25 GeV for a data subsample, 
corresponding to an integrated luminosity of 17 pb _1 , for which the track reconstruction 
efficiency at low momentum was smaller due to the operating conditions of the CTD [40]. 
For the untagged D° candidates, the incorrect assignment of the pion and kaon masses to 
the two tracks produces a wider reflected signal. This reflection was estimated from the 
D° candidates with a AM tag and normalised to the ratio of numbers of D° without and 
with AM tag; it was then subtracted from the untagged D° candidates. 

Figure [1] shows the M(Kir) distributions for untagged D° candidates after the reflection 
subtraction and for tagged D° candidates. The distributions were fitted simultaneously 
assuming that both have the same shape and are described by the "modified" Gaussian 
function [2]: 
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Gauss mod (M , a) oc e [-o.5^+V(i + o. 8 .-)^ (1) 

where x = \[M(K~tt + ) — Mo]/a\. The background shape in the fit was described [2] by 
the form [A + B ■ M(K-n + )] for M(K~n + ) > 1.86 GeV and [A + B ■ M(K-n + )} ■ exp{C • 
[M(K-tt + ) - 1.86]} for M{K~ti + ) < 1.86 GeV. The free parameters A, B and C were 
assumed to be independent for the two M(K~ir + ) distributions. Clear signals are seen at 
the nominal value of M(D°) [41]. The number of untagged (tagged) D° mesons yielded 
by the fit was N unt ^(D°) = 7996 ± 488 (N^^D ) = 1970 ± 78). 

5.2 Reconstruction of "additional" D* + mesons 

The D* + —> D %+ events with p T (D* + ) > 3 GeV and \r](D* + )\ < 1.6 can be considered as 
a sum of two subsamples: events with the D° having pt{D°) > 3 GeV and |?7(.D )| < 1.6, 
and events with the D° outside of that kinematic range. The former sample is a subset 
of D° mesons reconstructed with AM tag, as discussed above. The latter sample of 
additional D* + mesons was obtained using the same D° — > K~ir + decay channel. In each 
event, pairs of tracks with p? > 0.4 GeV were combined to form a D° candidate. Only 
combinations with invariant mass 1.80 < M(Kn) < 1.92 GeV were considered. The D° 
candidates were required to have pt(D°) < 3.0 GeV or \r](D°)\ > 1.6. A third track with 
Pt > 0.2 GeV, with charge opposite to that of the kaon track and assumed to have the 
pion mass, was combined with the D° candidate to form an additional D*^ candidate. 
Here again the cut value was px > 0.25 GeV for the data subsample for which the track 
reconstruction efficiency at low momentum was smaller. 

Figure [2] shows the distribution of the mass difference AM = M{K-ktx s ) — M(Ktt) for 
additional D* candidates. A clear signal is seen around the nominal value of M(D* ± ) — 
M(D°) [41]. The combinatorial background under the signal was estimated from the mass- 
difference distribution of the wrong-charge combinations, in which both tracks associated 
to the D° candidate have the same charge and the third track has opposite charge. The 
number of reconstructed additional D* + mesons was calculated by subtracting the wrong- 
charge AM distribution, after normalising it to the number of right-charge candidates in 
the region 0.150 < AM < 0.170 GeV. The subtraction was done in the signal region 
0.143 < AM < 0.148 GeV and yielded N add (D* + ) = 317 ± 26. 

5.3 Reconstruction of D + mesons 

The D + mesons were reconstructed from the decay channel D + — > K~7r + 7T + . The anal- 
ysis for this meson was done using the e + p data sample only, where the D + third level 
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trigger logic was implemented. In each event, two tracks with the same charge and a third 
track with opposite charge were combined to form a D + candidate. The tracks with the 
same charge were assigned the pion mass and required to have Pt( 7] ") > 0.5 GeV. For the 
remaining track, the kaon mass was assumed and Pt{K) > 0.7 GeV was required. Combi- 
natorial background was further suppressed by requiring cos8*(K) > —0.75. Background 
from D* + and — > <f)7r + with <p — > K + K~ was suppressed by applying suitable cuts [2]. 

Figure [3] shows the M(K~n + n + ) distribution for the D + candidates. A clear signal is seen 
at the nominal value of D + mass [41]. The mass distribution was fitted to a sum of the 
modified Gaussian function (Eq. ([]])) describing the signal and a linear function describing 
the non-resonant background. The number of reconstructed D + mesons yielded by the 
fit was N(D + ) = 4785 ± 501. 

5.4 Reconstruction of mesons 

The mesons were reconstructed from the decay channel — > cf)Ti + with <ft — > K + K~ . 
In each event, two tracks with opposite charges were each assigned the kaon mass and 
combined to form the candidates if the invariant mass M(KK) was within ±8 MeV 
of the nominal <fi mass [2]. A third track, assumed to be a pion, was combined with 
the (f) candidate, yielding the Df candidate. Only tracks with Pt{^) > 0.5 GeV and 
Pt{K) > 0.7 GeV were considered. To reduce the combinatorial background further, a 
cut was applied in the angle, 0*(tt), between the pion in the KKir rest frame and the 
KKtt line of flight in the laboratory frame, cos 9* (it) < 0.85, and in the angle, 0'(K), 
between one of the kaons and the pion in the KK rest frame, | cos 3 9'(K)\ > 0.1. 

Figure H] shows the M(K + K~n + ) distributions for the candidates after all cuts. A 
clear signal is seen at the nominal mass [41]. There is also a smaller signal around 
the nominal D + mass as expected from the decay D + — > </>7r + . The mass distribution 
was fitted to a sum of two modified Gaussian functions (Eq. ([T|)) describing the signals 
and an exponential function describing the non-resonant background. The number of 
reconstructed mesons yielded by the fit was N(D+) = 647 ± 80, for pt(D+) > 3 GeV 
and N(D+) = 773 ± 96, for p T (Df) > 2 GeV. 

6 Acceptance corrections 

The acceptances were calculated using the Rapgap 2.08 [42] MC model and checked with 
Herwig 6.3 [43]. The Rapgap MC model was interfaced with Heracles 4.6.1 [44] in 
order to incorporate first-order electroweak corrections. The generated events were then 
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passed through a full simulation of the detector, using Geant 3.13 [45], processed and 
selected with the same programs as used for the data. 

The MC models were used to produce charm by the BGF and the resolved photon pro- 
cesses. In the latter the virtual photon behaves as a hadron-like source of partons, one of 
which interacts with a parton of the initial proton. The CTEQ5L [46] and GRV-LO [47] 
PDFs were used for the proton and the photon, respectively. The charm-quark mass 
was set to 1.5 GeV. Both the Rapgap and Herwig MCs use LO matrix elements with 
leading-logarithmic parton showers. Charm fragmentation is implemented using either the 
Lund string fragmentation (in Rapgap) or a cluster fragmentation [48] model (in Her- 
wig). It was checked that both MC samples, Rapgap and Herwig, give a reasonable 
description of the data for DIS and D-meson variables when compared at detector- level. 

For a given observable Y, the production differential cross sections were determined using 

da_ _ N(D) 
dY ~ A-C-B-AY' 

where N(D) is the number of reconstructed D mesons in a bin of size Ay. The recon- 
struction acceptance A takes into account small admixtures in the reconstructed signals 
from other decay modes, migrations, efficiencies and QED radiative effects for that bin, C 
is the integrated luminosity and B is the branching ratio [41] for the decay channel used 
in the reconstruction (see Table [1]). The total production cross sections were determined 
using 

N(D) 
A-C-B" 

where N(D) and A are now for the whole kinematical range of the measurement. 

The reconstruction acceptances were calculated with Rapgap and vary depending on the 
particle and the kinematic region of the measurement. For 1.5 < Q 2 < 1000 GeV 2 , 0.02 < 
y < 0.7, transverse momenta pr(D , D + ) > 3 GeV, pt(D+) > 2 GeV and pseudorapidity 
\rj(D)\ < 1.6 the overall acceptances were ps 42%, rs 26% and « 17% for D°, D + , and D+ 
mesons, respectively. 

The relative 6-quark contributions, predicted by the MC simulation using branching ratios 
of 6-quark decays to the charmed hadrons measured at LEP [49,50], were subtracted from 
all measured cross sections. The subtraction of the 6-quark contribution reduced the 
measured cross sections by 3.1% for the D° and D + and 4.3% for the Df and changed 
the measured charm fragmentation ratios and fractions by less than 1%. 
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7 Charm- meson production cross sections 



Charm-meson cross sections were calculated using the reconstructed signals for the process 
ep -> eDX in the kinematic region 1.5 < Q 2 < 1000 GeV 2 , 0.02 < y < 0.7, p T {D) > 3 GeV 
(for the D+ also p T (Df) > 2 GeV) and \rj(D)\ < 1.6. 

The systematic uncertainties presented in this and the following sections will be discussed 
in Section [Til The third set of uncertainties quoted for the measured cross sections 
and charm fragmentation ratios and fractions are due to the propagation of the relevant 
branching-ratio uncertainties. 

The following cross sections were measured: 

• the production cross section for D° mesons not originating from the D* + — > D°ir^ 
decays, hereafter called untagged D° mesons, is: 

a nntag (D°) = 5.56 ± 0.35(stat.)+°;|(syst.) ± 0.10(br.) nb; 

• the production cross section for D° mesons originating from the D* + — > D°7r+ decays: 

o-tagpo) = i 78 ± 0.08(stat.)+°;^( sys t.) ± 0.03(br.) nb; 

• the production cross section for all D° mesons: 

a(D°) = a unUs (D°) + a ta,g (D°) = 7.34 ± 0.36(stat.)i°;^(syst.) ± 0.13(br.) nb; 

• the production cross section for additional D* + mesons: 

a add (D* + ) = 0.518 ± 0.046(stat.)^;°^(syst.) ± 0.01(br.) nb. 

The production cross section for D* + mesons in the kinematic range pt{D* + ) > 3 GeV 
and \v(D* + )\ < 1.6, a kin (D* + ), is given by the sum a* dd (D* + ) + a tag (D°) /B D *+^ D o n +: 

(7 kin (D*+) = 3.14 ± 0.12(stat.)^;^(syst.) ± 0.06(br.) nb; 

• the production cross section for D + mesons: 

a(D + ) = 2.80 ± 0.30(stat.)lau(syst.) ± 0.10(br.) nb; 
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• the production cross section for Df mesons with p T (D+) > 3 GeV: 



a(D+) = 1.27±0.16(stat.) 



+0.11 

-0.06 



(syst.) 



+0.19 
-0.15 



(br.) nb; 



• the production cross section for mesons with p T (D+) > 2 GeV: 



a 2 (D+) = 2.42 ± 0.30 (stat.) 



+0.30 
-0.14 



(syst.) 



+0.35 
-0.27 



(br.) nb. 



8 Charm fragmentation ratios 

In this section, the ratio of neutral to charged D-meson production rates, R u /d, the 
strangeness-suppression factor, j s , and the fraction of charged D mesons produced in 
a vector state, P^, are presented in the kinematic range 1.5 < Q 2 < 1000 GeV 2 , 0.02 < 
y < 0.7, pt(D) > 3 GeV and \rj(D)\ < 1.6. 

8.1 Ratio of neutral to charged D-meson production rates 

Neglecting influences from decays of heavier excited D mesons, the ratio of neutral to 
charged D-meson production rates is given by the ratio of the sum of D*° and direct D° 
production cross sections to the sum of D* + and direct D + production cross sections. 
This ratio can be calculated as [2] 



Using the measured cross sections, the ratio of neutral to charged D-meson production 
rates is 



The measured R u /d value agrees with unity, i.e. it is consistent with isospin invariance, 
which implies that u and d quarks are produced equally in charm fragmentation. 

In Tableland Fig. [5], this measurement is compared with those obtained in DIS by the HI 
collaboration [4], in photoproduction [2] and in e + e~ annihilations [51]. All measurements 
agree within experimental uncertainties. 




a untag (L>°) 



a(D+) + a ta g( J D )' 



Ru/d = 1.22 ± 0. 11 (stat. )±g;8i (syst.) ± 0.03(br.). 



10 



8.2 Strangeness-suppression factor 



The strangeness-suppression factor for charm mesons is given by the ratio of twice the 
production rate of charm-strange mesons to the production rate of non-strange charm 
mesons. All D* + and D*° decays produce either a D + or a D° meson, while all D* + 
decays produce a meson [41]. Thus, neglecting decays of heavier excited charm- 
strange mesons to non-strange charm mesons, the strangeness-suppression factor can be 
calculated as the ratio of twice the cross section to the sum of D° and D + cross 
sections. 

For the comparison of the inclusive _D + , D° and D* + cross sections with each other, the 
equivalent phase-space treatment [2] was used. The equivalent D + and D° cross sections 
were defined as the sum of their direct cross sections, i.e. D + and D° not coming from 
D*, and the contribution from D* + and D*° decays [2]: 



a^(D+) = a(D+) + a add (D* + )-(l-B D ^ D o n+ ), 

cr cq (D°) = a wata,g (D°) + a tag (D°) + a add (D* + ) ■ B D *+^ D o n+ + a add (/r ). 

It was also assumed that the production cross section for "additional" D*° mesons, pro- 
ducing D° mesons outside of the nominal kinematic range, is <j add (D*°) = R u /^a add (D* + ). 
The strangeness-suppression factor is then given by 



2a(Df) _ 2a(D s 

Is 



a^(D+) + a c i(D°) a(D+) + a xnskae (D°) + a ta s(D°) + cr* dd (D*+) ■ (1 + R u/d ) 
Using the measured cross sections, the strangeness-suppression factor is 



7s = 0.225 ± 0.030(stat.)t°;^ 8 7 (syst.)l[i;^(br.). 

Table [3] and Fig. [5] compare this measurement with the values measured in photoproduc- 
tion [2], in DIS by the HI collaboration [4] and in e + e~ annihilations [51]. All measure- 
ments agree within experimental uncertainties. The large branching-ratio uncertainties 
are dominated by the uncertainties of the Df — > (pir + branching ratio. 



8.3 Fraction of charged D mesons produced in a vector state 

Neglecting influences from decays of heavier excited D mesons, the fraction of charged D 
mesons produced in a vector state, P^, is given by the ratio of vector to (vector+pseudoscalar) 
charm meson production cross sections. The following relation holds [2]: 



11 



a(D+) + a ta g(D°) + a* dd (D*+) ' 

Using the measured cross sections, the fraction of charged D mesons produced in a vector 
state is 



P v d = 0.617 ± 0.038(stat.)±aoo9(syst.) ± 0.017(br.). 

The measured value is smaller than the naive spin-counting prediction of 0.75. Table H] 
and Fig. [5] compare this measurement with the values measured in photoproduction [2], 
in DIS by the HI collaboration [4] and in e + e~ annihilations [51]. All the measurements 
are consistent. 



9 Charm fragmentation fractions 

The fraction of c quarks hadronising as a particular charm meson, f(c — > D), is given by 
the ratio of the production cross section for the meson to the sum of the production cross 
sections for all charm ground states that decay weakly, cr gs . In addition to the measured 
D°, D + and charm ground states, the production cross sections of the A+ baryon 
and of the charm-strange baryons S+, S° and should be included in the sum. The 
production rates for the latter are expected to be much lower than that of the A+ due to 
strangeness suppression. 

The cr(A+) was estimated using the corresponding fragmentation fraction measured in 
e + e~ [51], f(c — > A+ ) e + e -, by the relation below: 

/(c^A c +) e+e -=a(A c +)Ar gs (2) 

The uncertainty of this procedure was estimated by using f(c—> A+) obtained in photo- 
production [2], and considering the uncertainty in f(c — > A+) e + e - [51]. 

The relative rates for the weakly-decaying charm-strange baryons were estimated from 
the non-charm sector following the LEP procedure [52]. The measured S~/A and Q~ /A 
relative rates are (6.65 ± 0.28)% and (0.42 ± 0.07)%, respectively [41]. Assuming equal 
production of S° and S~ states and that a similar suppression is applicable to the charm 
baryons, the total rate for the three charm-strange baryons relative to the A^ state is 
expected to be about 14%. Therefore, the estimated A+ production cross section was 
scaled by a factor 1.14 in the sum of the production cross sections. An error of ±0.05 was 
assigned to the scale factor when evaluating systematic uncertainties. 
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Using the equivalent D° and D + cross sections [2], a gs is given by 

a gs = a c \D + ) + cr cq (D°) + a(D+) + a(A+) • 1.14, 
which can be expressed as 

a gs = a(D + )+a mitas (D )+a tas (D°)+a add {D* + ) ■ (1 + R u/d ) + a(£>+) + <t(A+) -1.14. (3) 
Using the measured cross sections and combining Eqs. (j3J) and (J2J) yields 

a gs = 13.7 ± 0.6 (stat.)tj^ (syst.) ± 0.6 (br.) nb. 
The fragmentation fractions for the measured charm ground state are given by 

f(c^D+) = a c \D + )/a gs = [a(D + )+a* dd (D* + )-(l-B D * + ^ D o n+ )}/a gs , 

f(c^D°) = a^(D°)/a g8 

= [<t«***(D ) + a tas (D°) + a add (L>*+) • (R u/d + B D , + ^ + )]/a & , 

f{c^D+) = a(Df)/a gs . 

The fragmentation fraction for the D* + is 

/(c ^ D *+) = a kh \D* + )/a gs = [a^(D )/B D ^ D o w+ + cr axid (D* + )]/a gs . 

The open-charm fragmentation fractions, measured in the kinematic region 1.5 < Q 2 < 
1000 GeV 2 , 0.02 < y < 0.7, p T (D) > 3 GeV and \r](D)\ < 1.6, are summarised in Table E 
and Fig. The results are compared with the values obtained in photoproduction [2], in 
DIS by the HI collaboration [4] and in e + e~ annihilations [51]. All the measurements are 
consistent. A Monte Carlo study [2] suggested that the measured charm fragmentation 
ratios and fractions are close to those in the full Pt(D) and rj(D) phase space. 

The hadronisation fraction into untagged D°, needed in the next section for comparisons 
with theory, was: 



_>£>(,) = f( c ^ D °)- f ( c ^ D *+)B D . + -. D0 ^ 
= (a untas (D°) + <7 add (£>* ))/<7 flS 
= 0.450 ±0.020(stat.) lH(syst.) ^(br.), 

where the equivalent phase-space treatment and a add (-D* ) were considered as in Sec- 
tion El 
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10 Cross sections and pQCD comparisons 



For the cross sections presented in Section [7J the predictions from the HVQDIS program 
are a(D°) = 7.90 nb, a untag (D°) = 5.88 nb, (?{D + ) = 2.82 nb and a 2 (Df) = 2 AO nb, 
with uncertainties around 15%, dominated by the input PDF and the mass of the charm 
quark. They are in good agreement with the data. 

The differential cross sections for untagged D° (the D° mesons coming from D* + are 
already included in the previous ZEUS publication [3]), D + and as a function of Q 2 , 
x, Pt{D) and rj(D) are shown in Figs. 0, [7] and [8] and given in Tables O and [7J The cross 
sections in Q 2 and x both fall by about three orders of magnitude in the measured region. 
The cross-section da / dpr(D) falls by two orders of magnitude with increasing p T (D). The 
cross-section da/dr](D°) shows a mild increase with increasing rj(D°); for the D + and Df 
no statistically significant dependence with rj(D) is observed. 

Figures [6], [TJ and [8] show also the corresponding NLO calculations implemented in the 
HVQDIS program as well as their uncertainties (Section 111.21) . All the differential cross 
sections measured are well described by the NLO calculation. 

11 Systematic uncertainties 

11.1 Systematic uncertainties of measurements 

The systematic uncertainties of the measured cross sections and fragmentation ratios and 
fractions were determined by changing the analysis procedure and repeating all calcula- 
tions. 

In the measurement of fragmentation ratios and fractions the following groups of system- 
atic uncertainty sources were considered (Table [8]): 

• {S\} the model dependence of the acceptance corrections was estimated using the 
HERWIG MC sample; 

• {#2} the uncertainty of the beauty subtraction was determined by varying the fe-quark 
cross section by a factor of two in the reference MC sample; 

• {$3} the uncertainty of the tracking simulation was obtained by varying all momenta 
by ±0.3% which corresponds to the uncertainty in the magnetic field; and by changing 
the track momentum resolution and the angular resolution by ^0% of their values. 
The asymmetric resolution variations were used since the MC signals typically had 
somewhat narrower widths than observed in the data; 
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• {64} the uncertainty in the CAL energy scale was studied by varying in the MC the 
energy of the scattered e~ by ±1% and the energy of the hadronic system by ±3%; 

• {5$} the uncertainties related to the signal extraction procedures were studied as 
follows: 

— the cuts on the minimum pt for the it and K candidates were independently raised 
and lowered by 10% from their nominal values, 

— the cut on the minimum p? for the 7r s was raised and lowered by 0.02 GeV (for 
a tag (D°), a unta %(D°), a add (£>*+)), 

— the AM signal region was widened symmetrically by 0.003 GeV (for a tas (D°), 
a untag (D°), a add (D*+)), 

— the M(Kir) signal region was widened and narrowed symmetrically by 0.01 GeV 
(for a add (D*+)), 

— the wrong-charge background normalisation region was changed to 0.152 < AM < 
0.168 (for a add ( J D*+)); 

• {5q} the uncertainties of the luminosities of the e~p (±1.8%) and e + p (±2.25%) data 
samples were included, taking into account their correlations; 

• {£7} the uncertainty in the estimate of cr(A+) (see Section [9]); 

• {Sg} the uncertainty in the rate of the charm-strange baryons (see Section [9]); 

Contributions from the different systematic uncertainties were calculated and added in 
quadrature separately for positive and negative variations. Correlated systematic uncer- 
tainties largely cancel in the calculation of the fragmentation ratios and fractions. 

For the total and differential cross-section measurements discussed in Section [10] and 
those used for the extraction of F$° (Section fT2l) . further sources of systematics were 
studied [3,5,53], {S 9 }: 

• the cut on y e was changed to y e < 0.90; 

• the cut on yjB was changed to ?/jb > 0.03; 

• the cut on 5 was changed to 42 < 5 < 70 GeV; 

• the cut on |Z vcrtC x| was changed to |Z ver tex| < 45 cm; 

• the cut on E e > was changed to E e > > 11 GeV; 

• the excluded region for the impact position of the scattered electron in the RCAL was 
increased by 1 cm in each direction; 

• the electron method was used, except for cases when the scattered electron track was 
reconstructed by the CTD. In the latter case, the DA method, which has the best 
resolution at high Q 2 , was used. 
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These estimations were made in each bin in which the differential cross sections were 
measured. In addition, for the lowest x bin of the differential cross section of untagged 
D°, the systematic error accounted also for instabilities in the signal extraction, not en- 
countered in any other bin. The overall systematic uncertainty was determined by adding 
the individual uncertainties in quadrature. Typically 5$ was below 4%. The uncertainty 
on the luminosity measurement was not included in the systematic uncertainties for the 
differential cross sections. 



11.2 Uncertainties on theoretical predictions 

The NLO QCD predictions for D meson production are affected by the systematic uncer- 
tainties listed below. Typical values are quoted for the total cross section: 

• the ZEUS PDF uncertainties were propagated from the experimental uncertainties of 
the fitted data (±5%). As an alternative parametrisation in the FFNS, the CTEQ5F3 
PDF was used in HVQDIS with a charm mass of 1.3 GeV (—2%); 

• the charm mass was changed simultaneously in the PDF fit and in HVQDIS by 
=F0.15 GeV (lf%). The largest effect was at low p T {D); 

• the scale was changed to 2^JQ 2 + 4m;? an d to max(y / Q 2 /4 + mf, 2m c ) (_tg%); 

• the Jetset fragmentation as implemented in the previous analyses [3, 5] was used 
instead of the Peterson fragmentation (+5% to +20%). The largest deviations were 
observed for D° and D + particles at the lowest Q 2 and x. 



12 Extraction of F 2 CC 

The open charm contribution, F 2 CC , to the proton structure-function F 2 can be defined in 
terms of the inclusive double-differential cc cross section in x and Q 2 by 



d 2 a c5 (x, Q 2 ) 27ra 



2 



dxdQ 2 xQ A 



{[1 + (1 - y) 2 }Ff (x, Q 2 ) - y 2 F?( Xj Q 2 )}. 



In this paper, the cc cross section is obtained by measuring the untagged D°, D + and Df 
production cross sections and employing the measured hadronisation fractions /(c — > D). 
Since only a limited kinematic region in Pt{D) and r)(D) is accessible, a prescription for 
extrapolating to the full kinematic phase space is needed. 

As reported in Section [TjU the measured differential cross-sections are well described in 
the probed kinematic region. Therefore the following relation was used to extract F 2 CC : 



1(3 



tfte, Qi) = ai ' me i ep ^^ K,U^ Qi), (4) 

where crj imea s is the cross section in the bin i in the measured region of Pt{D) and rj(D) 
and (Ti it heo is the corresponding cross section evaluated with HVQDIS. The value of -F^theo 
was calculated in FFNS from the NLO coefficient functions [27] using the same values 
of parameters as in the calculation of cr^thco- The cross sections cx iimcas (ep — > DX) were 
measured in bins of Q 2 and y (Table [9]) and F 2 CC was quoted at representative Q 2 and x 
values for each bin (Table [TOl . The _F 2 CC measurements obtained from each D meson were 
combined into a single set of measurements; the result is also shown in Table [Ull 

The extrapolation factors from the measured Pt{D) and t](D) ranges to the full phase 
space, as estimated using HVQDIS, were between 17 at low Q 2 and 2.5 at high Q 2 for the 
D° and D + measurements. For the _D + , the lower pt requirement leads to lower extrapo- 
lation factors between 5.6 and 1.9. They are all shown in Table [TUl The uncertainty from 
the branching ratios was estimated by changing each branching ratio independently in the 
calculation by ±1 standard deviation and adding in quadrature the resulting variations 
ofFf (1^%). 

The following uncertainties of the extrapolation prescription of Eq. (j4j) have been evalu- 
ated: 

• using Jetset instead of the Peterson fragmentation yielded changes of ~ +28%, +15% 
and +5% for the data points at the lowest, middle and largest Q 2 ranges, respectively; 

• changing the charm mass by ±0.15 GeV consistently in the HVQDIS calculation and 
in the calculation of F£ c theo led to differences in the extrapolation of ±5% at low x, 
low Q 2 ; the value decreases rapidly at higher x and higher Q 2 ; 

• the uncertainty in the ZEUS NLO PDF fit led to uncertainties in the extracted values 
of _F 2 CC typically less than 1%; 

• the extrapolation factors were evaluated using the CTEQ5F3 proton PDF yielding 
differences compared to the nominal factors of ~ +10%, +6% and +1% for the lowest, 
middle and largest Q 2 ranges, respectively. 

The combined F£° measurements are shown in Fig. [9] The quadratic addition of the 
extrapolation uncertainties is shown independently as a band. Also shown in Fig. [9] is 
the previous measurement [3] and the ZEUS NLO QCD fit. The two sets of data are 
consistent!^. The prediction describes the data well for all Q 2 and x. The uncertainty 



The previous data were measured at Q 2 = 4, 18 and 130 GeV 2 and have been shifted to Q 2 = 4.2, 
20.4 and 111.8 GeV 2 using the ZEUS NLO QCD fit. 
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on the theoretical prediction is that from the PDF fit propagated from the experimental 
uncertainties of the fitted data. 

13 Summary and conclusions 

The production of the charm mesons D* + , D°, D + and D£ has been measured with 
the ZEUS detector in the kinematic range 1.5 < Q 2 < 1000 GeV 2 , 0.02 < y < 0.7, 
p T (D* + ,D°,D + ) > 3 GeV, p T (D+) > 2 GeV and \rj(D)\ < 1.6. 

The cross sections have been used to determine the charm fragmentation ratios and frac- 
tions. The ratio of neutral to charged D-meson production, R u /d, is compatible with 
unity, i.e. it is consistent with isospin invariance, which implies that u and d quarks are 
produced equally in charm fragmentation. The strangeness-suppression factor in charm 
fragmentation, j s , was measured to be about 20%. The fraction of charged D mesons 
produced in a vector state, P^, was found to be smaller than the naive spin-counting 
prediction of 0.75. The fraction of c quarks hadronising as D* + , D°, D + and -D+ mesons 
have been calculated. The measured R u /d, Is, P$ an d open charm fragmentation fractions 
are consistent with those obtained in charm photoproduction and in e + e~ annihilation. 
These measurements generally support the hypothesis that fragmentation proceeds inde- 
pendently of the hard sub-process. 

The measured D°, D + and differential cross sections were compared to the predictions 
of NLO QCD with the proton PDFs extracted from inclusive DIS data. A good description 
was found. 

The double-differential cross section in y and Q 2 has been used to extract the open 
charm contribution to F 2 , by using the NLO QCD calculation to extrapolate outside the 
measured Pt{D) and r](D) regions. The F% c values obtained from the different D mesons 
agree with previous results where a D* + meson was tagged. 
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Decay mode 


Branching ratio 


D° -> K~7T + 


0.0380 ±0.0007 


D* + -> D°7T+ 


0.6770 ± 0.0050 




0.0951 ±0.0034 


D+ -»■ -> K + K~TT + 


0.0216 ±0.0028 



Table 1: Branching ratios of the D-meson decay modes [41]. 
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Ru/d 


ZEUS (DIS) 


1.22 ± 0.11(stat.)iaoi(syst.) ± 0.03(br.) 


ZEUS ( 7 p) [2] 


1.100 ± 0.078(stat.)lH^(syst.)iHl9(br.) 


combined e + e~ data 


1.020 ± 0.069(stat. © syst.)l°;°^(br.) 


HI (DIS) [4] 


1.26 ± 0.20(stat.) ± 0.11(syst.) ± 0.04(br. © theory) 



Table 2: The ratio of neutral to charged D-meson production rates, R u /d- The 
e + e~ values are taken from [2]; they are an update of the compilation in [51] using 
the branching-ratio values of [54]- The measurements in this paper, ZEUS (DIS), 
took the values for all the branching ratios involved from [41], the rest of the quoted 
measurements took them from [54]. 





7, 


ZEUS (DIS) 


0.225 ± 0.030(stat.)+Ho 8 7 (syst.)lH26(br.) 


ZEUS ( 7 p) [2] 


0.257 ± 0.024(stat.)18:ow(syBt.)ig:8S(br.) 


ZEUS (tp) [11] 


0.27 ± 0.04(stat.)+°;^( sys t.) ± 0.07(br.) 


combined e + e~ data 


0.259 ± 0.023(stat. © syst.)l^g(br.) 


HI (DIS) [4] 


0.36 ± 0.10(stat.) ± 0.01(syst.) ± 0.08(br. © theory) 



Table 3: The strangeness-suppression factor in charm fragmentation, 7 S . The 
e + e~ values are taken from [2]; they are an update of the compilation in [51] using 
the branching-ratio values of [54]- The measurements in this paper, ZEUS (DIS), 
took the values for all the branching ratios involved from [41], the rest of the quoted 
measurements took them from [54]- 
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pd 

V 


ZEUS (DIS) 


0.617 ± 0.038(stat.)^aoo9(syst.) ± 0.017(br.) 


ZEUS ( 7 p) [2] 


0.566 ± 0.025(stat.)^;™I(syst.)i°;^(br.) 


combined e + e~ data 


0.614 ± 0.019(stat. © syst.)+°;°|(br.) 


HI (DIS) [4] 


0.693 ± 0.045(stat.) ± 0.004(syst.) ± 0.009(br. © theory) 



Table 4: The fraction of charged D mesons produced in a vector state, . The 
e + e~ values are taken from [2]; they are an update of the compilation in [51] using 
the branching-ratio values of [54]- The measurements in this paper, ZEUS (DIS), 
took the values for all the branching ratios involved from [41], the rest of the quoted 
measurements took them from [54]. 
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ZEUS [jp) [2\ 


uombined 








Pt{-L> ) > O bre V 


Pt\L > ) > o.o Lrev 


e + e~ data [51] 








\f]\U)\ < l.b 










stat. syst. br. 


stat. syst. 


stat.© syst. br. 


total 






0.216 ±0.019 


0.217 ±0.014 ±g;gj| 


0.226 ±0.010 to™ 


0.203 ±0.026 


f(c~ 


-> D°) 


0.605 ±0.020 l^l^ 


0.523 ± 0.021 ±g;gjf 


0.557 ±0.023 ±o:Sxs 


0.560 ±0.046 


f(c- 


-Dt) 


o.o92 ±0.011 ±85S±S:^ 


0.095 ± 0.008 ±g;gg| 


0.101 ±0.009 to°ol 4 


0.151 ±0.055 


f(c~ 


* D*+) 


0.229 ± 0.011 i^i^Z 


0.200 ± 0.009 to^Q 


0.238 ±0.007 ±g;gg| 


0.263 ±0.032 



Table 5: The fractions of c quarks hadronising as a particular charm hadron, 
f(c — > D). The fractions are shown for the D + , D° and charm ground states 
and for the D* + state. The measurements in this paper, ZEUS (DIS), took the 
values for all the branching ratios involved from [41], the rest of the quoted mea- 
surements took them from [54]. 



untagged D D_ 



c^r Din vjev j 


An 1 AC) 2 
CIO 1 UA^l 


A 

^stat 
^IID/ VortiV ) 


A 

^syst 


An /ah 2 

CIO / CLLJ 


A 

^stat 
^IlD/ *ort!V ) 


A 

^syst 


1 5 5 


56 


+0 08 


+0.09 
-0.03 


275 


±0 066 

_I_W . WUU 


+0.056 
-0.044 


5, 15 


0.141 


±0.008 


+0.011 
-0.006 


0.089 


±0.014 


+0.003 
-0.010 


15, 40 


0.044 


±0.005 


+0.003 
—0.002 


0.016 


±0.003 


+0.002 
—0.001 


40, 1000 


0.0012 


±0.0002 


+0.0002 
-0.0001 


0.0007 


±0.0002 


+0.0001 
-0.0001 




LIU j LLJb 


^-*stat 
l nD J 


A 

"syst 


LIU I LLJb 


A 

^stat 
l nD J 


A 

"syst 


000091 0001 


1 51 97 


±2543 


+1137 
-5789 


4470 


±1 572 

_l L (J 1 


+1078 
-1112 


0.0001, 0.0005 


4162 


±304 


+457 
-231 


2295 


±327 


+157 
-199 


0.0005, 0.001 


1476 


±195 


+128 
—72 


1049 


±178 


+39 
— 125 


0.001, 0.1 


20.7 


±2.3 


+2.7 
-1.0 


10.5 


±2.2 


+1.4 
-0.5 


Pt\U ) Din ijrev ) 


An l Arn i r>o\ 
cio I cipx[J-J ) 


A 

^stat 


A 

^syst 


An 1 Ar> l n+\ 

cio l cipTyiy ) 


A 

^stat 

fnh 1 CA\J\ 
^IIU/ VjrtjV ) 


A 

^syst 


o.u, o.o 


O.lO 


QQ 


+0.28 
-0.12 


l.Ul 


-4-0 AA 


+0.45 
-0.21 


3.5, 4.5 


1.93 


±0.20 


+0.12 
-0.11 


0.77 


±0.14 


+0.09 
-0.06 


4.5, 6.0 


0.78 


±0.11 


+0.05 
-0.08 


0.49 


±0.08 


+0.05 
-0.02 


6.0, 20. 


0.051 


±0.009 


+0.004 
-0.003 


0.028 


±0.007 


+0.002 
-0.001 


t](D) bin 


do/d V {D°) 


Astat 

(nb) 


^syst 


do/d7]{D + ) 


Astat 

(nb) 


^syst 


-1.6, -0.6 


1.18 


±0.19 


+0.13 
-0.10 


0.65 


±0.11 


+0.08 
-0.09 


-0.6, 0.0 


1.59 


±0.19 


+0.10 
-0.11 


1.20 


±0.25 


+0.15 
-0.22 


0.0, 0.6 


2.05 


±0.22 


+0.18 
-0.14 


1.06 


±0.21 


+0.08 
-0.15 


0.6, 1.6 


2.31 


±0.37 


+0.09 
-0.20 


0.74 


±0.23 


+0.22 
-0.07 



Table 6: Measured differential cross sections for D° not coming from a D* + 
(left), and D + (right) as a function of Q 2 , x, Pt{D) and t](D) for 1.5 < Q 2 < 
1000 GeV 2 , 0.02 < y < 0.7, p T (D) > 3 GeV and \q{D)\ < 1.6. The estimated 
b-quark contribution of 3.1 % has been subtracted. The statistical and systematic 
uncertainties are shown separately. The D° (D + ) cross sections have a further 1.8% 
(3.6%) uncertainty from the D° — > K~n + (D + — > K~ix + ix + ) branching ratios. 
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Q 2 bin ( GeV 2 ) 


da/dQ~ 


A 

^stat 


A 

^syst 








(nb/ GeV ) 




1.5, 


rr 
O 


0.31 


±0.07 


+0.08 
-0.05 


5, 


15 


U.Uyz 


±U.Ulo 


+0.004 
-0.017 


10, 




0.016 


±0.005 


+0.004 
-0.003 




1 nnn 


00025 


±0 nnm n 

_I_W .UuulU 


+0.00008 
-0.00004 


a; 


bin 


7 / J 

da /ax 


A 

^stat 

(nb) 


A 

^syst 


0.000021, 


0.0001 


4982 


±1967 


+1354 
-1333 


0.0001, 


0.0005 


Z (00 


i A An 

±44o 


+65 
-644 


0.0000, 


0.001 


934 


±250 


+118 

-155 


0.001, 


0.1 


0.1 


1 1 r- 
±1.0 


+0.8 
-0.6 




bin ( GeV) 


da/dp T (Dl) 


A 

^stat 

(nb/ GeV) 


A 

^syst 


2.0, 


2.5 


1.00 


i n rr) 

±U.0z 


+0.36 
-0.50 


2.5, 


3.0 


0.62 


±0.22 


+0.14 
-0.11 


3.0, 


3.5 


0.59 


i n m 

±0.21 


+0.08 
-0.12 


3.5, 


4.5 


0.00 


i r\ i i 
±U.ll 


+0.05 
-0.05 


4.5, 


6.0 


0.20 


±0.05 


+0.02 
-0.01 


6.0, 


20. 


0.011 


±0.004 


+0.002 
-0.001 


rj(Df) bin 


da/d V (Df) 


^stat 


A SyS t 








(nb) 




-1.6, 


-0.6 


0.94 


±0.24 


+0.11 
-0.26 


-0.6, 


0.0 


0.57 


±0.15 


+0.14 
-0.04 


0.0, 


0.6 


0.76 


±0.18 


+0.06 
-0.09 


0.6, 


1.6 


0.85 


±0.22 


+0.17 
-0.12 



Table 7: Measured differential cross sections as a function of Q 2 , x, pt{D^) 
and r](D+) for 1.5 < Q 2 < 1000 GeV 2 , 0.02 < y < 0.7, p T (D+) > 2 GeV and 
\r](Df)\ < 1.6. The estimated b-quark contribution of 4-3 % has been subtracted. 
The statistical and systematic uncertainties are shown separately. The cross sec- 
tions have a further 13% uncertainty from the — > <f)7r + — > K + K~ir + branching 
ratio. 
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°1 ~~ °9 


°1 ~~ °8 


A, 
°1 


°2 


°3 


°4 


fir 

°5 


°6 


°7 


°8 






y/o) 


y/o) 




y/o) 


y/o) 


y/o) 


y/o) 


y/o) 


y/o) 


(7 untaS (£> ) 


+5.8 


+5.4 


+2.5 


+1.7 


+1.3 


+1.5 


+3.5 


+2.2 






-4.7 


-4.1 


-0.0 


-3.3 


-0.9 


-0.6 


-0.4 


-2.1 






a{D+) 


+6.6 


+3.2 


+0.6 


+1.5 


+1.0 


+1.4 


+0.0 


+2.3 






-5.0 


-4.6 


-0.0 


-3.1 


-1.9 


-0.6 


-1.5 


-2.2 








+9.1 


+8.8 


+0.0 


+0.0 


+1.2 


+0.5 


+8.0 


+2.2 






-7.4 


-7.2 


-2.0 


-4.0 


-0.0 


-0.4 


-5.1 


-2.1 






0- tag (Z>°) 




+5.6 


+0.0 


+1.8 


+2.3 


+1.5 


+4.0 


+2.2 








-4.6 


-1.5 


-3.5 


-1.1 


-0.3 


-0.8 


-2.1 






a add (D*+) 




+9.9 


+0.0 


+1.9 


+4.2 


+2.3 


+8.2 


+2.2 








-8.9 


-2.9 


-3.8 


-0.2 


-0.5 


-7.2 


-2.1 






a kin (D* + ) 




+5.7 


+0.0 


+1.8 


+2.5 


+1.6 


+4.0 


+2.2 








-4.7 


-1.8 


-3.6 


-0.8 


-0.3 


-1.1 


-2.1 










+8.9 


+2.8 


+2.2 


+4.0 


+0. 1 


+6.8 


+2.2 








-4.9 


-0.0 


-4.5 


-0.0 


-0.1 


-0.0 


-2.1 






Ru/d 




+4.3 


+2.7 


+0.0 


+1.4 


+0.1 


+3.0 


+0.1 








-1.4 


-0.0 


-0.1 


-1.4 


-0.1 


-0.3 


-0.1 






Is 




+7.9 


+1.8 


+0.5 


+4.1 


+0.6 


+6.4 


+0.0 








-3.0 


-0.0 


-1.2 


-0.4 


-1.3 


-2.4 


-0.0 






pi 




+2.8 


+0.0 


+0.1 


+1.4 


+0.2 


+2.4 


+0.0 








-1.4 


-1.3 


-0.3 


-0.0 


-0.0 


-0.0 


-0.1 






er gs 




+10.3 


+1.2 


+1.7 


+1.2 


+1.4 


+2.4 


+2.2 


+9.3 


+0.4 




-4.3 


-0.0 


-3.4 


-0.5 


-0.5 


-0.5 


-2.1 


-0.9 


-0.4 


/(c-D+) 




+1.1 


+0.0 


+0.3 


+0.1 


+0.1 


+0.2 


+0.1 


+0.9 


+0.4 




-9.2 


-0.7 


-0.2 


-1.7 


-0.2 


-2.9 


-0.1 


-8.5 


-0.4 






+1.1 


+0.0 


+0.0 


+0.5 


+0.2 


+1.8 


+0.0 


+0.9 


+0.4 




-8.6 


-0.0 


-0.0 


-0.4 


-0.1 


-0.0 


-0.0 


-8.5 


-0.4 


yuntag ( c £)0) 




+2.0 
-8.6 


+1.0 
-0.0 


+0.1 
-0.0 


+0.5 
-0.9 


+0.2 
-0.1 


+1.4 
-0.2 


+0.0 
-0.0 


+0.9 
-8.5 


+0.4 
-0.4 






+7.1 


+1.6 


+0.5 


+3.7 


+0.5 


+5.7 


+0.0 


+0.9 


+0.4 




-8.9 


-0.0 


-1.0 


-0.3 


-1.2 


-2.1 


-0.0 


-8.5 


-0.4 


/(c^D*+) 




+2.8 


+0.0 


+0.1 


+1.5 


+0.2 


+2.2 


+0.0 


+0.9 


+0.1 




-9.1 


-2.9 


-0.1 


-0.3 


-0.0 


-1.5 


-0.0 


-8.5 


-0.4 



Table 8: The systematic uncertainties resulting from 5i-5g, from 61-6$, and from 
5\ to 5s independently (see text) for the charm hadron cross sections and charm 
fragmentation ratios and fractions. 
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untagged D 



Q 2 bin ( GeV 2 ) 


y bin 


a A stat A syst (nb) 


1.5, 9 


0.18, 0.70 
0.02, 0.18 


1.50 ±0.19 t° it 
1.03 ±0.19 toil 


9, 44 


0.20, 0.70 
0.02, 0.20 


1.02 ±0.14 t° f 9 
1.02 ±0.14 ^11 


44, 1000 


0.02, 0.70 


1.03 ±0.19 to.H 


D + 


Q 2 bin ( GeV 2 ) 


y bin 


& A stat A syst (nb) 


1.5, 9 


0.18, 0.70 
0.02, 0.18 


0.63 ±0.14 t°of 
0.65 ±0.13 f° f 8 


9, 44 


0.20, 0.70 
0.02, 0.20 


0.52 ±0.11 t f 
0.44 ±0.11 t° ° t 


44, 1000 


0.02, 0.70 


0.61 ±0.25 to° l 


Df 


Q 2 bin ( GeV 2 ) 


y bin 


a A stat A syst (nb) 


1.5, 9 


0.18, 0.70 
0.02, 0.18 


0.97 ±0.23 toil 
0.56 ±0.15 t£og 


9, 44 


0.20, 0.70 
0.02, 0.20 


0.61 ±0.15 toit 
0.29 ±0.07 jftJl 


44, 1000 


0.02, 0.70 


0.20 ±0.11 t° ° l 



Table 9: Measured cross sections for D not coming from a D* + , D + and 
in each of the Q 2 and y bins for p T (D°,D + ) > 3 GeV, pr(Df) > 2 GeV and 
1 77 I < 1.6. The estimated b-quark contribution of 3.1 % for D° and D + and 4-3% 
for Df has been subtracted. The statistical and systematic uncertainties are shown 
separately. The D°, D + and Df cross sections have further 1.8%, 3.6% and 13% 
uncertainties from the D° — > K~ir + , D + — > K^tt + ti + and Df — > (pir + — > K + K~tt + 
branching ratios, respectively. 
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untagged D 



Q 2 (GeV 2 ) 


X 


-^2 A s tat A sys t 


A cx trap factor 


4.2 


0.00013 
0.00061 


0.141 ±0.017 to° ll 
0.090 ±0.017 to° u 


+0.048 o n 
-0.013 °- y 

+0.036 i 7 
-0.006 1 ' 


20.4 


0.00062 
0.00281 


0.320 ±0.044 t° f 2 l 
0.156 ±0.021 igloos 


+0.061 4 q 
-0.020 

+0.041 r o 
-0.004 °-° 


111.8 


0.00676 


0.217 ±0.039 lg;gf? 


+0.014 9 r 
-0.002 z - d 



Q 2 (GeV 2 ) 


X 


Astat A syst 


A cx t r ap factor 


4.2 


0.00013 
0.00061 


0.123 ±0.025 1°;°^ 
0.109 ±0.020 Iffi 


+0.037 o q 
-0.011 °- J 
+0.039 i 7 
-0.007 1 ' 


20.4 


0.00062 
0.00281 


0.331 ±0.067 toTj 
0.130 ±0.039 tof l 


+0.066 4 q 

-0.021 

+0.030 r o 
-0.003 °-° 


111.8 


0.00676 


0.293 ±0.124 tftJS 


+0.021 9 r 
-0.003 



Q 2 (GeV 2 ) 


X 


FT A sta t A syst 


Aextrap factor 


4.2 


0.00013 
0.00061 


0.221 ±0.044 to^ll 
0.075 ±0.017 ±8:°n 


+0.036 j o 
-0.016 

+0.019 r fi 
-0.004 °' u 


20.4 


0.00062 
0.00281 


0.470 ±0.100 t&i™ 
0.100 ±0.022 ±2:oot 


+0.037 9 o 
-0.017 

+0.013 2 q 

-o.ooi 


111.8 


0.00676 


0.179 ±0.058 ±g:o3B 


+0.013 i q 
-o.ooi 



Combined 



Q 2 (GeV 2 ) 


X 


■^2° A s t a t A syst 


A ex trap 


4.2 


0.00013 
0.00061 


0.144 ±0.014 t° f 2 5 

0.090 ±0.010 ±S:SSS 


+0.045 
-0.013 
+0.029 
-0.005 


20.4 


0.00062 
0.00281 


0.341 ±0.035 t °oit 
0.132 ±0.014 t f t 


+0.063 
-0.021 
+0.024 
-0.001 


111.8 


0.00676 


0.211 ±0.032 


+0.013 
-0.002 



Table 10: The extracted values of from the production cross sections of D° not 
coming from D* + , D + and Df and the combination of them at eachQ 2 andx value. 
The statistical, systematic and extrapolation uncertainties are shown separately. 
The values of the extrapolation factor used to correct the full Pt{D) andrj(D) phase 
space are also shown. All the extracted F| c values have a further ±2.7% —4.1% 
uncertainty from the D° — > K~-k + , D + — > A'~7r + 7r + and Df — > <p7i + — > K + K~n + 
branching ratios and the f(c—> A+) value. 
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Figure 1: The M(K~tt + ) distributions (dots) for (a) the D° candidates without 
AM tag, obtained after reflection subtraction, and for (b) the D° candidates with 
AM tag. The lowest and highest mass bins are affected by the trigger selection. 
The solid curves represent a fit to the sum of a modified Gaussian function and a 
background function. 
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Figure 2: The distribution of the mass difference, AM = M(K~7r + 7rf) — 
M(K~-k + ) , for the "additional" D* + candidates (dots). The histogram shows the 
AM distribution for wrong-charge combinations. For illustration, the solid curve 
represents a fit to the sum of a modified Gaussian function and a background func- 
tion. 



32 



ZEUS 



> 

CD 



CD 
CL 

CO 

a 
o 

H — ' 

CC 
c 

E 
o 
O 



16000 



14000 



12000 



10000 



8000 



6000 



4000 



2000 




1.7 



— 


— 




- 
- 






• ZEUS (65 pb" 1 ) 


* 


Gauss mod + Backgr. 




1.5<Q 2 <1000 GeV 2 




0.02<y<0.7 




P T (D + )>3 GeV 




h(D + )|<1.6 




N(D + ) =4785 ±501 





1.8 



1.9 



2 2.1 
M(Ktot) (GeV) 



Figure 3: The M(K~ir + 7r + ) distribution for the D + candidates (dots). The 
lowest and highest mass bins are affected by the trigger selection. The solid curve 
represents a fit to the sum of a modified Gaussian function and a linear background 
function. 
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Figure 4: The M(K + K~tt + ) distribution for (a) candidates with pr{D+) > 3 
GeV and (b) Df candidates with pt{D^) > 2 GeV. The lowest mass bins are 
affected by the trigger selection. The solid curves represent fits to the sum of two 
modified Gaussian functions and an exponential background function. The first 
peak in both distributions is from D + decaying through the same channel. 
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Figure 5: (a) The ratio of neutral to charged D-meson production rates, R u /d, 
the strangeness-suppression factor in charm fragmentation, 7 S; and the fraction of 
charged D mesons produced in a vector state, P^. (b) The fractions of c quarks 
hadronising as D + , D° and Df charm ground-state mesons, as D* + mesons and 
as A+ baryons. The inner error bars show the statistical uncertainties and the 
outer bars show the statistical and systematic uncertainties added in quadrature. 
The measurements have further uncertainties coming from the different branching 
ratios involved; their magnitudes are shown in Tables [|| and [7] for R u /d, 7s and 
Py, respectively, and in Tabled for the fractions. 
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Figure 6: Differential cross sections for D° not coming from D* + as a function 
of Q 2 , x, p T (D°) and rj(D°) compared to the NLO QCD calculation of HVQDIS. 
The inner error bars show the statistical uncertainties and the outer bars show the 
statistical and systematic uncertainties added in quadrature. The lower and upper 
NLO QCD predictions show the estimated theoretical uncertainty of the HVQDIS 
calculations. The data points have a further 1.8% uncertainty from the D° — > 
K~7i + branching ratio. 
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Figure 7: Differential D + cross sections as a function ofQ 2 , x, px(D + ) andrj(D + ) 
compared to the NLO QCD calculation of HVQDIS. The inner error bars show 
the statistical uncertainties and the outer bars show the statistical and systematic 
uncertainties added in quadrature. The lower and upper NLO QCD predictions 
show the estimated theoretical uncertainty of the HVQDIS calculations. The data 
points have a further 3.6% uncertainty from the D + — > K~ir + 7r + branching ratio. 
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Figure 8: Differential Df cross sections as a function ofQ 2 , x, pt(D^) andrj(Df) 
compared to the NLO QCD calculation of HVQDIS. The inner error bars show 
the statistical uncertainties and the outer bars show the statistical and systematic 
uncertainties added in quadrature. The lower and upper NLO QCD predictions 
show the estimated theoretical uncertainty of the HVQDIS calculations. The data 
points have a further 13% uncertainty from the — > <p7i + — > K + K~n + branching 
ratio. 
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Figure 9: The measured F| c as a function of x for three Q 2 bins. The cur- 
rent data are compared with the previous ZEUS measurement [3]. The data are 
shown with statistical uncertainties (inner bars) and statistical and systematic un- 
certainties added in quadrature (outer bars). All measured _F 2 CC values have a further 
+2.7% -4-1% uncertainty coming from the current experimental uncertainty from 
the D° — > K~ir + , D + — > K~7r + 7r + and Df — > <pir + — > K + K~ir + branching ratios 
and the f(c—> A+) value. The shaded band corresponds to the estimated theoret- 
ical uncertainty in the extrapolation. The lower and upper curves show the ZEUS 
NLO QCD fit [27] uncertainty propagated from the experimental uncertainties of 
the fitted data. 
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